Trinomial Tree Construction 


A trinomial tree based method is presented for pricing exotic options. The model is based on a 
combination of techniques. that is, a tree generation technique and an appropriate backward 


induction pricing technique. 


Since the volatility parameter in the SDE is of a piecewise constant form, the tree generation 
techniques may, in some cases, construct trees that are non- recombining. In the worst case, then, 
the space complexity of the tree generation techniques is proportional to the exponential of the 


number of time slices in the tree. 


Let 0=1f, <:--<t, =T bea partition of the time interval [0,7]. Furthermore suppose that the 


underlying security follows piecewise geometric Brownian motion, in the sense described below, 


over the interval [0,7]. Specifically, assume that the underlying security can be modeled as a 


process, {s (t)|t €[0, T}} , which, under the risk neutral probability measure, satisfies a stochastic 


differential equation (SDE) of the form 


dS(t) = ut) S(t)dt + o(t)S(t)dwWit), t €[0,7], (1) 


where {w(n)|e €[0, T}} is standard Brownian motion. Here s(t) and o(t) are deterministic 


functions of the piecewise constant form 


u, t €[0,¢,), o,, t €[0,t,), 
L(t) = and o(t) = 
Ly> belt aby | On> b-eltyasly 


Each method includes a technique for constructing, based on the SDE (1), an appropriate tree of 
discrete prices of the underlying security. Each such technique uses a mathematical result, 
described below, for ensuring that branching probabilities from each tree node are appropriate 


(i.e., probabilities, for each node, must be non-negative and sum to one). 


Consider a tree node, @, ata time slice, t;, where 0<i< N; furthermore, assume that the 
logarithm of the price of the underlying security at this node is equal to log Sold. We assume 


that node @ branches into three nodes, at time slice f, 


i+1? 


with respective logarithm of the price 
of the underlying security of the form (A+ l)AlogSnew, (A)AlogSnew, and 
(A-l)AlogSnew where 4 €Rand Alog Snew > 0. 


Here (A)A log Snew is the value that, among all tree nodes at time ¢,,,, is closest to 


i+1? 
2; 
i+] 


A k (oy 
log Sold + u,,,At;., where Hiwt = Hist — 2 


and At,,, =¢,,, —¢,; furthermore, 


(A +lAlog Snew and (A —1)Alog Snew are values for the two nodes closest to the node with 
value (A)AlogSnew. Next we associate with node @ a discrete random variable, Y, which 


takes the values 


(A + l)Alog Snew, with probability p,, 
Y =) A(Alog Snew), with probability p,,, 
(A —l)Alog Snew, with probability p, 


We seek to determine p,, p,, and p,, above, so that the mean and variance of the discrete 


random variable Y match those of the continuous random variable 


log Sold + f,,,At,,, +O, 


ay . (obtained by solving the SDE (1), with initial condition 


log S(t, ) = log Sold, for the time interval [t,,¢,,,]). 


By matching mean and variances as described above, and by ensuring that the probabilities sum 


to one, we obtain the following system of linear equations 


Dp, (A+ DA log Snew + p,,A(A log Snew) + p, (A — lA log Snew = log Sold 
* Pat 5 ? 


p, (A +1) (Alog Snew)* + p,, (A)? (A log Snew)? + 
p,(A-1)’ (Alog Snew)’ =o,,,’ At 


i+l i+1 


+ (logSold + ji,,,At.)?, 


Pu +D,, + Dy = 1, 


for the unknowns p,,p, andp,,. By algebraic manipulation, the linear system of equations 


above is equivalent to 


Py —Pa =B-A, 
p, (1+ 2A) + p, (1-2) = A’ +B-7V, (3) 


Py + Pm + Pa =1, 


Where 


NT (4 ) 
= a 
(Alog Snew)* 


2 O74; 


and 


_ log Sold + 14,,,At,,, (4b) 
7 A log Snew 


Notice that while the system of equations above has a unique solution, we have no guarantee that 
P,>Pm and p, will be non-negative. Next we determine a condition on A to ensure that 


PusPmoPa 20), 


Recall that the branching rule from node @ implies that 


1 ‘ 1 
(A+ aa log Snew = log Sold + w,,,At,,, =(A- aA log Snew. (5) 


Dividing (5) by Alog Snew and substituting B for the right hand side of (4b), we have 


1 1 
aA ts >BreA- >? which we can rewrite as 


5 (6) 


for some x €[0,1]. Solving (3) and substituting the right hand side of (4b) for B, we obtain 


Ae 
oe ae any 
Py a gtx— A, (7) 
Ae ae 
Pax 5 + 5 = 


where x €[0,1]. Notice that the right hand side of (7) has no dependency on 1. 


Analyzing the right hand side of (7) over the range x €[0,1], we obtain the condition 


(8) 


on A, which ensures that p,,),,,P 4 29. Notice that (8) yields an equivalent condition on 


A log Snew , that is, 


(9) 


20;, V At,, 
=—G < Alog Snew < 20,,,,/At,,; - 


To summarize, for an arbitrary tree node on an arbitrary time slice, appropriate branching 


probabilities are given as the solution of (2) provided that condition (9) holds. In Appendix B 


we examine the sensitivity of solutions to (2) with respect to perturbations in the values of 


AlogSnew and J. 


In this section we present the techniques for generating a tree appropriate for pricing the barrier 


options described in Section 2. We consider single barrier options first. 


Suppose that we have constructed a tree up to time t,, where 1<i<N. To expand the tree to 


the next time slice, we first define, at time ¢,,,, an appropriate partition for the logarithm of the 


underlying security; then, using this partition, we determine the children and associated 


probabilities of all nodes at time f;. 


Note that, by an appropriate partition for the logarithm of the underlying security at time f,,,, 


we mean a partition such that the inter-node spacing is equal to AlogSnew where Alog Snew is 


chosen (as in Section 3.1.1) so that branching probabilities are non-negative. Next we describe 
how to construct such a partition. Then we discuss how to determine the branching and 


corresponding probabilities for nodes at time f,. 


To define a partition at time ¢,,, (with uniform spacing, Alog Snew , which satisfies the 
inequality (9)), we first determine whether, for some nodes at time ¢; , there is a branch that 
crosses the barrier at time ¢,,,. This determination is made by checking certain conditions, 


defined in Appendix A, based on the branching rule (for nodes on the old time slice to nodes on 


the new time slice) defined. I 


f we determine that the barrier will be crossed at time t, 


i+1? 


we generate a partition by placing a 
node on the actual barrier (i.e., either H/,(¢,,,) or H,(t,,,)) and all other nodes offset from the 
barrier by integer multiples of AlogSnew. Otherwise, if the barrier will not be crossed, we 
define an artificial barrier at time f,,, (see Appendix A); we then generate a partition by placing 


a node exactly on the artificial barrier and all other nodes offset by integer multiples of 


Alog Snew from this barrier (here the artificial barrier simply acts as a point of reference for 
generating the partition). We use for Alog Snew a value close to the upper bound, 20,,,,/At;,, 


in the inequality (9). 


Once an appropriate partition has been defined at time f,,,, we then determine, according to the 


i+]? 
branching rule presented in Section 3.1.1, the children of each node at time ¢,. Suppose that a 
particular node at time ¢, (with value for the logarithm of the underlying security equal to 


log Sold ) branches to nodes at time f¢,,, with values for the logarithm of the underlying security 


i+ 


equal to a,, a,, and a,, respectively, where a, >a,, >a, . Then, from (2), appropriate 


u? 


branching probabilities are determined by solving the system of linear equations 


Dj, Gy + Pin Gm + Pada = log Sold + f1,,,At 


i+]? 


Pa, + P,Q, + pa; = Gig Abs + (log Sold + f,,,At,,,)° 
| ite Ute ot oF el 


for the unknowns p,,p, and p,,. The numerical conditioning of the system of linear 


equations above should be checked, however, to ensure the accuracy of the computed solution. 


The tree construction technique for double barrier options is based on a similar approach as for 


single barrier options. That is, if the tree has been constructed up to time ¢;, an appropriate 


partition for the underlying security is defined at time ¢,,,; then branches and associated 


probabilities are determined for nodes on the old time slice. We describe these techniques next. 


Suppose that the tree has been generated up to time ¢,, where 1<i< JN. If neither barrier is 


crossed at time f, 


i+]? 


then an artificial barrier is defined, at time ¢,,,, and used (as a point of 


reference) to generate an appropriate partition. Otherwise, a partition is defined that places 
nodes on both barriers (see below). Once a partition is defined, branching and corresponding 


probabilities are determined as above. 
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